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Abstract 

Breaking of isotropy and Lorentz boost invariance in the dynamics of second-gen- 
eration leptons would lead to direction-dependent changes in the lifetimes of charged 
pions. This would make the intensity of a neutrino beam produced via pion decay a 
function of the beam orientation. The experimental signature of this phenomenon — 
sidereal variations in the event rate at a downstream neutrino detector — has already 
been studied, in searches for Lorentz-violating neutrino oscillations. Existing analyses of 
MINOS near detector data can be used to constrain the flavor- diagonal Lorentz violation 
coefficients affecting muon neutrino speeds at roughly the 10~^ level. 



baltschuOpliysics . sc . edu 



1 Introduction 



One exotic topic in elementary particle physics that has received a significant amount 
of attention of late is the possibility of that the fundamental laws of nature may not 
be perfectly invariant under Lorentz and CPT symmetries p]. Thus far, there is no 
compelling experimental reason to doubt that these symmetries are exact. However, there 
are important theoretical reasons to be interested in the topic. Many prospective theories 
of quantum gravity may involve the breaking of these symmetries, at least in certain 
regimes. The mechanisms for the symmetry breaking include spontaneous breaking in 
string theory 0, E] and elsewhere [3], effects in loop quantum gravity |S1 E] and non- 
commutative geometry [TJ [8] , Lorentz violation through spacetime- varying couplings [9l 
[To] , and anomalous breaking of Lorentz and CPT symmetries pT] . 

The discovery of any kind of Lorentz or CPT violation would be extraordinarily 
important — a sure sign of new physics, with a fundamentally different structure from any- 
thing that has so far been observed. Such a discovery would provide a remarkable window 
onto a fundamentally new and different physical regime. This fact makes the subject area 
quite interesting, for both theoretical and experimental research programs. Most recent 
theoretical work is now performed within the context of effective quantum field theory. 
There is an extensively-studied effective field theory known as the standard model exten- 
sion (SME) that contains all possible translation-invariant but Lorentz- violating operators 
that may be constructed out of known standard model fields. 

In the SME, the effects of Lorentz-violating operators are parameterized by small 
tensor- valued coefficients fi2\ [13] . In scenarios in which the Lorentz violation arises spon- 
taneously, these tensor coefficients represent the nonzero vacuum expectation values of 
tensor-valued fields. Because CPT violation in a quantum field theory that is both stable 
and unitary automatically entails the existence of Lorentz violation as well [H] , the SME 
is automatically the general effective field theory for describing CPT-violating physics. 

As a systematic effective field theory, the SME contains, in principle, an infinite num- 
ber of possible Lorentz-violating operators. However, all but a small number of these 
operators are expected to be suppressed at modest experimental energies. For low-energy 
experiments, the most relevant subset of the general theory is the minimal SME, which 
involves only those forms of Lorentz and CPT violation that are gauge invariant, local, 
and superficially renormalizable. This minimal theory has become the standard frame- 
work used for parameterizing experimental Lorentz tests. Recent searches for Lorentz 
violation have included studies of matter-antimatter asymmetries for trapped charged 
particles [151 [ISl [IZj and bound state systems [ISl [12], measurements of muon prop- 
erties [20|, [21], analyses of the behavior of spin-polarized matter p2], frequency stan- 
dard comparisons [231 1211 |25l [26], Michelson-Morley experiments with cryogenic res- 
onators [271 [281 [291 [301 [31] , Doppler effect measurements [321 [33], measurements of neutral 
meson oscillations [331 IS51 [SSI [23 [3S1 131] , polarization measurements on the light from 
cosmological sources [3D1 El 1321 133] , high-energy astrophysical tests [331 [3S1 1301 133 I3E] , 
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precision tests of gravity |l9l |50], and others. The results of these experiments set con- 
straints on the various SME coefficients, and up-to-date information about most of these 
constraints may be found in |51] . 

The neutrino sector of the standard model is one of the most difficult to study, since 
neutrinos only interact weakly. However, there have been a number of analyses of how 
Lorentz violation could affect neutrino oscillation phenomena. The general framework for 
studying SME-driven neutrino oscillations was given in ^2]. Based on this work, there 
have been experimental analyses of the results from the MINOS [531 ISU [551 [56] , LSND [57] . 
IceCube [2H], Double Chooz [52], and MiniBooNE [BD] experiments. This paper will look 
at a different kind of neutrino Lorentz violation, one not involving oscillations; however, 
the ultimate physical observable is the same as in some of the oscillation analyses. This 
means that elements drawn from previous data analyses can immediately be carried over 
into the present work. 

Measuring flavor-diagonal Lorentz violation for neutrinos can be a challenge. The 
MINOS analyses were intended as a search for off-diagonal SME coefficients — for exam- 
ple, looking for oscillations — )■ along a baseline whose direction is changing slowly 
with the Earth's rotation. Particle oscillation experiments can be extremely sensitive to 
small effects, because they measure the quantum mechanically coherent build-up of tiny 
effects. In contrast, the conventional way to identify neutrino Lorentz violation that does 
not involve flavor-changing operators is to look for deviations from the expected neutrino 
velocity = 1 [HH E21 E31 El]- Measuring the velocities of neutral particles such as 
neutrinos is much more difficult than equivalent measurements with charged species; for 
charged leptons, it is possible to measure their velocities indirectly, through their inter- 
actions with photons, because it is the four-velocity that governs how a charge interacts 
with the electromagnetic field. 

This paper will continue to focus on those forms of Lorentz violation that can modify 
the neutrino velocities most pronouncedly. However, rather than direct measurements of 
Vu, the focus will be on other effects generated by the same SME coefficients. In particular, 
changes to the limiting neutrino velocity will also affect the rate at which neutrinos are 
produced in weak decays. 

This paper is organized as follows. Section [2] will discuss how Lorentz violation in the 
lepton sector affects the decay rates for charged pions. Section [3] will discuss the influence 
of the speciflc experimental geometry of the NuMI beam on pion decay observables. New 
bounds based on previous analyses of the MINOS data are presented in section |H followed 
by a discussion of the outlook for future experimental improvements. 

2 Lorentz Violation in Pion Decay 

Most of the SME coefficients for the first-generation fermions that make up everyday 
matter have been measured quite precisely. However, forms of Lorentz violation involving 
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only second-generation particles are quite a bit more difficult to constrain. Experiments 
that are sensitive to second-generation Lorentz violation must involve either unstable 
charged particles or weakly coupled neutrinos. Nonetheless, there are observables that 
depend rather sensitively on the SME coefficients for the second-generation species. In 
particular, there are flavor-diagonal effects that could produce the same kind of signature 
in the MINOS near detector as Lorentz-violating oscillations between and another 
species u^. These phenomena cause the initial intensity of the neutrino beam to depend 
on the direction in which the beam is pointing. Periodic changes in the neutrino beam 
intensity would then change the number of muon charged current events recorded in 
the the detector, giving the same experimental signature as Lorentz-violating neutrino 
oscillations. 

The most prominent reason that the NuMI beam intensity might vary with the sidereal 
orientation is that Lorentz violation for the second-generation leptons affects the lifetimes 
of the pions that decay to produce the beam. In any decay, the dispersion relation for the 
daughter particles has a strong effect on the decay rate, since it determines the amount 
of phase space available for the decay products. Modiflcations of the muon and neutrino 
energy-momentum relations will therefore change the rate at which pions decay to produce 
those particles. There are also changes, of comparable size, to the weak interaction matrix 
element for the process. 

The details of the decay rate have been worked out [65j , to leading order in the relevant 
Lorentz violation coefficients. The SME Lagrange density relevant for the leptons in this 
scenario is 

C = iL{Y + cl^-i^)D^L + iR{-i^' + Cj^-i^)D^R + mass terms, (1) 
where L and R are the left- and right-chiral lepton multiplets, 

, R=[^r]. (2) 

and Z)^ is the covariant derivative, containing the electromagnetic and weak gauge inter- 
actions. There are separate coefficients for the left-chiral and right-chiral leptons; these 
can be alternately expressed in terms of Lorentz violation coefficients d^^ and dy^, where 
= d^^ + (F^ and c"^ = d^^ — (F^. However, since only left-chiral particles undergo 
charged current weak interactions, only cl affects the decay rate of charged pions. The 
effects of flavor- diagonal minimal SME coefficients other than the c and d are suppressed 
at relativistic energies. 

Moreover, it is important to note that SU{2)l gauge invariance requires that the 
Lorentz violation coefficients {clY^ be the same for the left-chiral charged lepton and its 
corresponding neutrino. This means that it is not possible to have neutrinos with Lorentz- 
violating dispersion relations without having the same kinds of modiflcations also applying 
to the charged leptons — unless the chiral gauge invariance (which is responsible for such 
features as lepton universality) is abandoned. 
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The cl coefficients together form a traceless background tensor. The isotropic element 
of this tensor, (cl)oo, endows spacetime with a preferred frame structure. This leads to an 
isotropic violation of Lorentz boost symmetry. The other eight independent coefficients 
also break spatial isotropy as well as boost invariance. Although Lorentz violation is often 
accompanied by CPT violation, the effects of the cl coefficients are all CPT-even. 

The pion decay rate T depends, at leading order, on the isotropic coefficient (cl)oo, 
evaluated in the rest frame of the decaying pion. The dependence is [65] 




where Tq is the rest frame decay rate in the absence of Lorentz violation. The effects of 
(cl)oo are enhanced by the similarity in the pion and muon masses; using the physical 
values of these quantities, the decay rate is F = Fq [1 + 9.4(ci)oo]- 

In the pion rest frame, the total decay rate depends only on the (ci)oo in that frame. 
However, because decaying pions are typically moving relativistically relative to the lab- 
oratory, {cl)oo depends on the direction of that motion. This gives the pion lifetime 
observed in the laboratory a sidereal dependence, as the pion beam direction moves with 
the rotation of the Earth. Notably, the sidereal dependence of the decay rate (and thus 
the neutrino beam strength) is equally sensitive to violations of isotropy and violations of 
boost invariance. This is a sharp contrast with the sensitivities seen in many low-energy 
laboratory experiments, which are typically much more sensitive to anisotropy than boost 
violation. The reason for the parity in sensitivity in this case is that the decaying pions are 
moving relativistically in the laboratory, and this provides strong sensitivity to relativistic 
effects; as the Earthbound laboratory rotates, this provides a way of testing whether the 
behavior of the pions is the same when they are highly boosted in different directions. 

Experimental constraints on the SME coefficients are conventionally expressed using 
a set of Sun-centered celestial equatorial coordinate |66] . The origin of these coordinates 
lies at the center of the Sun. The Z-axis points parallel to the Earth's rotation axis; 
the X-axis points toward the vernal equinox point on the celestial sphere; and the Y- 
axis is determined by the right hand rule. Time in these coordinates is T; however, it 
is generally advantageous to choose a local time coordinate so that at = 0, the 
laboratory ^/-direction coincides with the F-direction in the Sun-centered frame. 

The isotropic {cl)oo in the pion rest frame is 



in terms of the pion velocity tv, its Lorentz factor 7,^, and the SME coefficients in the 
Sun-centered frame. The velocities associated with the movement of the Earthbound 
laboratory are all nonrelativistic, so is essentially the the pion's velocity relative to the 
lab. 

The linear combination {cl)tt + {cl){tj){vtt)j + (cl) jx(^^7r) j(^'7r)_ft: is the key quantity 
that determines the energy-momentum relation for ultrarelativistic leptons moving in the 





(3) 



(cl)oo = ll {cl)tt + {cl){tj){v^)j + {cl).jk{v^).j{vt,)k 



(4) 
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direction -0^. The dispersion relation for the highly relativistic decay products (as observed 
in the lab frame) is E = [1 - {cl)tt - {cl){t.j){v^)j - {cl)jk{v^)j{v.„)k\\p\- This takes 
into account relativistic beaming, which causes the decay products to be emitted along 
essentially the same direction as the parent pion. 

The dependence of F on the pion beam direction is not the only manifestation of 
anisotropy. In the pion rest frame, the decay is actually anisotropic. Decay products are 
more likely to emerge in some directions than others, but the increased rates of decays 
into certain directions compensate for decreased rates in other directions, leaving the total 
decay rate dependent only on (cl)oo (at lowest order). However, this kind of anisotropy is 
difficult to measure, because relativistic effects ensure that the decay products are beamed 
into a narrow pencil of angles around the direction of the pion's motion. It is worth 
noting that beaming is also responsible for the fact that it is only the dispersion relation 
for particles moving essentially in the {),r-direction that determines the laboratory frame 
decay rate. That this should be the case can be seen from the fact that if the decay rate 
were (cumbersomely) evaluated in the lab frame, the phase space available for the decay 
products would only depend on their energy-momentum relations within their narrow 
range of allowed angles around v.,,. 

While a full understanding of the decay anisotropy in the pion rest frame would be 
fairly interesting, it is beyond the scope of the present analysis, since this information 
is not required in order to place order of magnitude bounds on the flavor- diagonal cl 
parameters. In fact, calculating the angular dependence of the decay would require a 
substantially more elaborate generalization of the calculation in [65]. There are no ap- 
proximations apparent that could reliably be used to simplify such a calculation. Phase 
space estimates are sometimes useful in evaluations of the rates for Lorentz-violating pro- 
cesses. These estimates involve ignoring the effects of Lorentz violation on the dynamical 
matrix elements involved and studying merely how the Lorentz violation coefficients affect 
the phase space kinematics. This strategy is typically useful when looking at processes 
that are normally forbidden by Lorentz symmetry — that is, those for which the phase 
space available for the process is conventionally zero [Ill[67ll68]. It is not surprising that 
in such situations, the changes to the kinematics should predominate. However, the cal- 
culations in [65J show that for pion decay, the effects of the cl coefficients on the dynamics 
and kinematics are comparable in size; a phase space estimate of the pion decay rate will 
not be quantitatively accurate. 

Another simplification would apply if both the daughter particles were moving ultra- 
relativistically in the center-of-mass frame, as is the case for the vr"*" — )■ e"*" + Ve decay 
mode. In such a scenario, there is no dynamical or kinematical distinction between the 
two outgoing particles in the absence of Lorentz violation [because of the Lorentz- invariant 
theory's SU{2)l and CP symmetries]. The equivalence of the daughter particles, which 
travel in antipodal directions in the center-of-mass frame, would preclude any anisotropy 
that had different rates of neutrino production in the p- and — p-directions. This would 
rule out all odd multipole patterns; in particular, there could be no dependence on the 
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dipolar (cL)(oj) coefficients, and any anisotropy would have to come from the {cL)jk- Un- 
fortunately, for the decay to a muon-neutrino pair, the muon is not energetic enough for 
this simplification to occur. 

3 Effects of Experimental Geometry 

In a neutrino oscillation experiment like MINOS, a detector intercepts the neutrino beam 
near the beam's point of origin. If Lorentz violation in the neutrino sector leads to 
sidereal oscillations in the strength of the neutrino beam, this will be detectable using the 
near detector. In this section, we will focus specifically on understanding how to make 
such a measurement using the NuMI neutrino beam and the MINOS detector located at 
Fermilab, because there is more information available about possible sidereal oscillations 
for MINOS than for other high-energy neutrino experiments. 

The most elementary constraint based on the variability of the pion decay rate arises 
from the fact that there can be a threshold boost beyond which the decay 7r+ — )■ /i^ + 
becomes energetically impossible. This effect is nonperturbative in the cl coefficients, 
and so it was not included in the perturbative analysis that led to ([3]). Nonetheless, 
if > ^{ml - ml)/2[-{cL)TT - {.CL)(Tj)iv-K)j - {cL)jKiy-K)j{v^)K\, the greater-than- 
normal growth of the muon and neutrino energies as functions of \p \ means that there is 
not enough energy to produce daughter particles with the required momentum [69| ITO]. 
The primary pion decay channel becomes disallowed. 

Naturally, the existence of an upper threshold for the pion decay process could lead 
to sidereal variations in the beam. The nonoccurence of the decay above a certain pion 
energy will affect the neutrino intensity. As the beam direction rotates with the Earth, 
the threshold will vary, and a variable fraction of the expected population of high energy 
neutrinos will fail to be produced. However, a constraint derived from an analysis of this 
phenomenon appears to be less sensitive than one that can be extracted from an analysis 
of the decay rate based on the perturbative formula Q. 

The energy distribution of the neutrinos in the NuMI beam is somewhat complicated. 
However, in order to get an order of magnitude constraint on the coefficients for 
second-generation leptons, we shall treat the problem as the decay of a ~ 6 GeV pion, 
without regard to the direction of the decay neutrino. This energy is chosen to simulate 
a neutrino beam with a peak at ~ 3 GeV. A full analysis of the beam energy structure 
would introduce a number of additional complications. The degree of forward beaming 
depends on the progenitor pion energy, and this interacts nontrivially with the effects of 
Lorentz violation as well. While the total decay rate F depends only on the isotropic 
(cl)oo in the pion rest frame, the amplitudes for specific angular decay channels depend 
on the anisotropic coefficients contracted with the momentum vectors of the daughter 
particles. Certain neutrino production directions will be favored. Especially at lower pion 
energies, where beaming is less pronounced, this can affect the energy composition of the 
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NuMI beam. In contrast, at high energies, the decay neutrinos are strongly ahgned into 
the observable beam, and the sensitivity to Lorentz violation is simultaneously higher, 
because of a larger factor 7,^ appearing in the boosted expression for {cl)oo. 

Other experimental complications may be dealt with more explicitly. The result ([3]), 
in conjunction with the formula for (cl)oo, gives the instantaneous decay rate for pions 
moving in a particular direction. However, what is measured with an observation of the 
neutrino beam strength is not the decay rate but the total number of pions that manage 
to decay during the neutrino beam production process. This means that the length of the 
pion decay pipe affects the magnitude of any variations in the neutrino beam strength. If 
the pion decay rate, as measured in the laboratory frame, is 7,rr; then for pions moving 
at speed v-,, ^ 1 along a decay pipe of length D, the fraction that decay over the length of 
the pipe is P{D) = 1 — e^'^'"^^. The fractional difference in P{D), if there is a fractional 
change in the decay rate AF/Fo, is 

^AF^ 7Jo^ ... 




r=rn 



Note that for D 3> (77rFo)~^, essentially all the pions have time to decay, so a change in 
the decay rate does not affect the total beam strength. In contrast, if D ^ (7^Fo)~^, the 
decaying fraction is P{D) ^ '-j^J^D, and the fractional change in the number of decays 
is simply AF/Fq. As D increases, the sensitivity of the beam intensity is suppressed; 
however, there is also naturally an improvement in statistics at the detector, as the total 
number of the neutrinos in the beam grows. The optimal sensitivity for a monoenergetic 
beam will occur in the vicinity of D ~ (77rro)^- For the NuMI beam, the decay pipe 
length is 677 m long, chosen to be comparable to {;^t,Vq)^^ = 7^^(2.6 x 10"^ s) for pions 
with GeV energies. For a 6.0-GeV pion, the suppression factor 'y-^ToD/^e'^''^^^ — 1) is 
0.31. 

Finally, the most important influence that the specific geometry of the MINOS ex- 
periment has on the possibility of sidereal oscillations in the detector event rate comes 
from the details of how the pion decay direction changes with time. Let the detector be 
located at colatitude x ix = 42.18° for the MINOS near detector at Fermilab), and let 
[6, (f)) be spherical coordinates describing the angle between the beam direction and the 
local zenith direction {6 = 93.27°) and the azimuthal angle in the plane of the Earth's 
surface, measured starting eastward from south (0 = 203.91°). Then the beam direction 
at a local time = is [7T] 

= NiX + N^Y + N^Z (6) 
= (cos X sin 6 cos (p + sin x cos 6)X + (sin 6 sin (j))Y 

+ (— sinxsin6'cos0 + cosxcos^^)Z (7) 

= -0.715X- 0.405y + 0.571Z. (8) 

As the Earth rotates with sidereal frequency w®, the key quantity (cl)oo varies according 
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to 

(cl)oo = In [A + cos(w®T©) + Buj sin(w®Te) + A20J cos(2w©Te) + B2uj sin(2u;©Te)] 
where the coefficients are 



(9) 



Ao = icL)TT + N3{cL)iTZ)+Nl{cL)zZ + ^il-Nl)[{cL)xX + {cL)YY] (10) 
= Ni{cL)iTX)+NiN3{cL)(XZ)+N2{cL)iTY)+N2N3{cL)(YZ) (H) 

Bu, = -N2{cl)(tx)-N2N:,{cl)(^xz) + Ni{cl)(ty)+NiN^{cl)(yz) (12) 

= \{Nl-Nl){cL)^+N,N2{cL)iXY) (13) 

B2. = -iViiV2(cL)- + ^(iV2-iV|)(ci)(xy), (14) 

with (cl)_ = {cl)xx - {cl)yy- 

The ultimate experimental observable — the number of muon events recorded in the 
detector — will therefore be modulated at the frequencies and 2uj^. Collecting all the 
factors to determine the sensitivity, the fractional modulation in the normalized event 
rate will have amplitudes equal to the A and B coefficients times a sensitivity factor 
S composed of the Lorentz factor squared [from (jl])], the decay pipe factor, and the 
enhancement factor from ([3]): 

S = (7^ { ^ff ] { ^1 = 5.8 X 10^ (15) 

The largest part of the enhancement obviously derives from the highly relativistic nature 
of the progenitor pions. 



4 Conclusions From MINOS Data 

In order to constrain the coefficients for the second-generation leptons using MINOS 
data, a analysis of the observed beam intensity versus sidereal time is required. Fortu- 
nately, the core of the necessary analysis has already been performed by the experimenters; 
this was the basis of the papers [531 ES] , the conservative spirit of whose analyses will be 
followed here. The observables for the scenario with a direction-dependent rate of oscilla- 
tions z/^ — )■ Vx and the scenario with a similarly direction-dependent pion decay rate are 
the same: a change in the number of muon charged current events that indicate interac- 
tions of the beam neutrinos with the detector. Since [53| [55] found no significant evidence 
for sidereal oscillations in the number of muon neutrinos in the beam, there can be no 
evidence of a variation in the charged pion decay rate in these data sets either. Instead, 
there will be constraints on the anisotropic parameters. 
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Table 1: Bounds on the magnitudes of the individual cl coefficients for second-generation 
leptons, assuming only one coefficient at a time is nonvanishing. 



The MINOS analyses found no evidence of oscillations that reached the 3a level of 
significance. In fact, the observed signal levels were well below such a level. The analysis 
in [53] determined the Fast Fourier Transform power present in each quadrature mode at 
frequencies and 2u^, as well as additional higher harmonics. The level of statistical 
noise was characteristic of a Icr dispersion in these power values of 1.8 x 10~^, when the 
data was normalized to the total event rate (which was set by the rate of protons on 
target in the collisions that originate the pions). Thus a 3a bound on the normalized 
amplitudes SA and SB of the oscillations at the two lowest frequencies is 5.6 x 10~^. 
This corresponds to constraints on the separate A and B coefficients of approximately 
9.7 X 10"^ 

The results of [55] are similar but used an antineutrino data set with slightly less 
statistical power. Note that even in the presence of nonzero d coefficients, the antineutrino 
data provide complementary constraints on the same cl coefficients, which govern the 
dispersion relations for left-handed neutrinos and their right-handed antiparticles. 

Of the nine Lorentz- violating components of {ciY^, six can be constrained by looking 
at sidereal variations in the event rates at multiple neutrino experiments. The ones 
that cannot be constrained this way include {cl)tt^ which is purely isotropic, as well as 
(cl)(tz) and {cl)zzi which describe the effects of a preferred direction lying parallel to 
the Earth's axis. In the spherical harmonic notation of [72], the coefficients affecting 
only a single generation of leptons are equivalent to the {c^i)jm, which run up to highest 
weights up to j = 2; the six subject to constraint according to this method are those 
with m 7^ in celestial equatorial coordinates. For a single particular experiment, four 
linear combinations of the remaining six parameters may be constrained. These linear 
combinations are represented by the A and B coefficients. 

If the absence of sidereal oscillations in the signal is not due to a coincidental cancel- 
lation among the coefficients, it is possible to set order of magnitude constraints on the 
six separate \{cl){tx)\, |(cL)(ry)|, \{cl)-\, \{cl){xy)\, \{cl){xz)\ and \{cl){yz)\ coefficients. 
This may be done by setting all but one of the coefficients to zero, then finding the value 
of the remaining coefficient that is required to raise one of the ^ or ;B to an observable 
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3a level. The results are given in table [TJ 

These order of magnitude constraints on the anisotropy coefficients, at the 10~^ level, 
are comparable to bounds on time-independent deviations of v^, from 1 based on terrestrial 
experiments with muon neutrinos. The sidereally averaged neutrino speed measured in 
such experiments depends the six constrained here and also on the remaining coefficients 
{cl)tt, {cl){tz), and {cl)zz, as outlined in [65j. Each such speed measurement depends 
on a particular linear combination of the cl coefficients involved. 

In future searches for anisotropic neutrino oscillations, based on muon neutrino dis- 
appearance, it should be possible to add additional sensitivity to the flavor- diagonal co- 
efficients discussed here. The sensitivity to the cl coefficients that affect solely second- 
generation leptons derives from the way the form of Lorentz violation involved affect the 
decay of the pion progenitors of the beam neutrinos. The sensitivities of future mea- 
surements obviously depend most on the energies of the pions involved; greater 7^ values 
will generally improve constraints. The position and orientation of the beam have lesser 
influences, as does the length of the pion decay region. A very long decay pipe actually 
wipes out sensitivity in the lower-energy component of the neutrino beam; the beam in- 
tensity is independent of at a given energy if essentially all the pions in the relevant 
energy range have sufficient time to decay. It is possible that this kind of shift toward 
higher energies may actually be advantageous; although the signal is lost at lower energies, 
the sensitivity to the cl coefficients increases with 7^, and the quality of the relativistic 
beaming approximation also increases with energy. A full understanding of these effects 
will require more detailed understanding of the parent pion and daughter neutrino beams' 
energy structures. 

The current level of sensitivity is much lower than the sensitivity to SME parameters 
that induce oscillations. It also poorer than the ~ 10~^^ sensitivity that has been obtained 
from the observed absence of the photon decay reaction 7 — )■ /i~ -|- /i"*" for cosmic ray 
photons [73]; however, because the photon is a spin-1 state, the daughter particles in 
such a reaction would necessarily have parallel spins, which makes such constraints only 
sensitive to the spin- averaged c coefficients, rather than the left-chiral cl- On the other 
hand, the sensitivity of the new method discussed in this paper is comparable to that 
achieved for flavor-diagonal effects using direct time-of-flight measurements, and this new 
strategy can represent an important tool for constraining boost and isotropy violation in 
the difficult-to-access neutrino sector. 
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